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The neutron capture cross sections (n, y) of Bromine was performed with the time-of-flight (TOF) 
technique at the Back-n facility of the China Spallation Neutron Source (CSNS). Prompt y-rays 
originating from neutron-induced capture events were detected by four CeDe detectors. The pulse 
height weighting technique (PHWT) and the double-bunch unfolding method based on the Bayesian 
theory were used in the data analysis. To obtain reliable measurement results, background deduc- 
tions, normalization, and corrections were carefully considered. The multilevel R-matrix Bayesian 
code SAMMY was used to extract the resonance parameters in the resolved resonance region (RRR). 
The average cross sections in unresolved resonance region the was obtained from 10 to 400 keV. The 
experimental results were compared with the data of several evaluated libraries and previous ex- 
periments in the resolved and unresolved resonance region. The TALYS code was used to describe 
the average cross sections in the unresolved resonance region. Astrophysical Maxwell average cross 
sections (MACSs) of “8! Br from kT = 5 to 100 keV was calculated over a sufficiently wide range of 
neutron energies. At a thermal energy of kT = 30 keV, the MACS value for “Br 638464 mb is in 
good agreement with the KADoNiS recommended value. By contrast, the value of 293429 mb for 
8l Br is substantially higher than those of the evaluated database and the KADoNiS recommended 
value. 
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Fig. 1. Nucleosynthesis processes occurring in the region 
among Br. The s- and r-processes are indicated by black 
and blue arrows, respectively. 


and could be connected with the corresponding stellar 
objects and sites|1, 2]. The main s-process occurs in 
the He-rich inner shell of thermally pulsing asymptotic 
giant branch (AGB) stars and produces predominantly 
nuclei with mass numbers A>90. The weak component, 
which is responsible for the production of nuclei between 
iron and yttrium (56<A<90), takes place during convec- 
tive core He burning stage in massive stars|1]. Since the 
neutron exposure is small, the weak s-process flow can 
not overcome the bottleneck at the closed neutron shell 
N= 50[2]. 

The s-process components of the solar abundances 
in the Br-Kr-Rb region are characterized by the super- 
position of abundance contributions from the main s- 
process associated with thermally pulsing low mass AGB 
stars and from the weak s-process[3, 4]. However, the 
two components contribute in completely different ways. 
While the relative strength of the weak component de- 
creases significantly with mass number, that of the main 
component is rapidly growing. The complexity of this 
situation is further enhanced if one recalls that the weak 
and the main s-process exhibit two very different neu- 
tron capture regimes. For this reason, the branches in 
this regime are sketched in Fig. 1, represent the most 
problematic part of the entire s-process path. Quantita- 
tive investigations of these aspects have to rely on accu- 
rate stellar (n, y) cross section. Furthermore, lanthanum 
bromide detector are used in nuclear experiments for 
neutron and y-ray detection, accurate and reliable ex- 
perimental data on the neutron-induced reaction are re- 
quired for detectors design and optimization. 

As far as the available data for Br are concerned, pre- 
vious TOF experiments and works based on the acti- 
vation technique were mainly located in unresolved res- 
onance region (URR)[5, 6, 15]. Fig. 2 shows the pre- 
vious experimentally (n,y) cross sections of Bromine. 
M. Ohkubo et al(1981)[13] and R. L. Macklin (1988)[14] 
measured neutron resonance parameters of “Br and 
®1Br with electron accelerator neutron source. Our 
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Fig. 2. Previous experimentally measured cross sections of 
nat Br (nyy) compiled in the EXFOR database[7, 12, 15, 16]. 


work was performed at the back-streaming white neutron 
beam line (Back-n) of China Spallation Neutron Source 
(CSNS)[17-19, 21-23]. The deexcited y-rays are de- 
tected by four hydrogen-free deuterated benzene (C6De) 
liquid scintillator detectors[18, 24]. A °LiF—Si detector 
array is used as a neutron flux monitor to provide in- 
line neutron monitoring and the neutron flux profile[25]. 
The experimental conditions of the Back-n facility used 
in this study and the corresponding experimental setup 
are presented in Sec. II. The data analysis is provided in 
Sec. III, including the PHWT, double-bunch unfolding 
method, background deductions, normalization, and cor- 
rections. The theoretical description of our experimental 
results is discussed in Sec. IV. Finally, the conclusions 
are drawn in Sec. V. 


Il. MEASUREMENTS 


Measurements were made with the CeDe detector sys- 
tem in the Back-n facility of CSNS[28]. Neutrons are 
produced by slamming 1.6 GeV/c? proton beam, with 
double bunches per pulse, onto a tungsten target with a 
typical repetition rate of 25 Hz. The pulse width of each 
bunch is 41 ns, and the interval between the two bunches 
is 410 ns[29]. Along the neutron beam-line, there are 
two experimental stations: the near station (ES#1), 56 
m from the spallation target to the sample, and the far 
station (ES#2) with near 76 m neutron flight path[28). 
The detector system for the (n,y) reaction measurement 
consists of four CgD¢ detectors, the aluminum detector 
brackets, and an aluminum sample holder, as shown in 
Fig. 3. The CDe detectors are placed upstream of the 
sample, and the detector axis is at an angle of 110° rela- 
tive to the neutron beam[24]. A ŝLiF-Si detector array 
with 360 ug/cm? LiF neutron conversion layer and eight 
separated Si detectors was used for neutron flux moni- 
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Fig. 3. Photograph of the CgD¢ detector system in the Back- 
n facility of CSNS; a gold sample is placed in the aluminum 
sample holder. 


toring. Signals from these detectors are processed by a 
generalized full waveform digital data acquisition system, 
in which the flash analog-to-digital converters are based 
on the folding-ADC and FPGA techniques. Each chan- 
nel has a digital resolution of 12 bits and a sampling rate 
of 1 GSPS, corresponding to a time step of 1 ns/sample. 
The digital waveform data of the detectors are filtrated 
by a full-digital trigger system and then transferred to 
the CSNS computation center for long-term storage[25]. 


Our work was carried out at experimental station 
ES#2, the shutter and collimators had an inner diameter 
of 650+ 615+ 640 mm, resulting in a circular Gaussian- 
shaped beam profile with a diameter of around 40 mm 
at the sample position[29]. A thin foil of a cadmium 
absorber was placed at the front of the neutron shut- 
ter to absorb neutrons with an energy lower than 0.5 
eV for avoiding the overlap between consecutive neutron 
pulses. In addition, a Ta—Ag—Co filter with a total 
thickness of 1.0+0.4+1.4 mm was used to determine the 
in-beam y-ray background by employing the black res- 
onance method. A total of 5 samples were used in the 
measurements: (i) the KBr crystal; (ii) a 1°’ Au sample 
for experimental setup verification and flux normaliza- 
tion; (iii) a natural carbon sample and (iv) a lead sample 
to determine the backgrounds due to scattered neutrons 
and the in-beam y rays, (v) empty target target to de- 
termine sample-independent background. In the experi- 
ments, the samples were fixed on the aluminum sample 
holder of the CgDg detector system. The characteris- 
tics of the samples and the experimental duration are 
summarized in Table 1. 
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Ill. DATA ANALYSIS 


A. PHWT 


The efficiency of CgD¢ detectors for detecting the neu- 
tron capture events is depending on the deexcitation 
paths of the compound nucleus, which are too complex to 
be calculated[30]. The pulse height weighting technique 
(PHWT) is required in the measurements to manipulate 
response function of CgDg detector for y-rays. In a such 
way, detection efficiency (ey) for y-ray with energy of 
Ey will satisfy 


(1) 


where c is a proportionality coefficient. When ey is suf- 
ficiently low (ey < 1), the efficiency for detecting a cap- 
ture event (€e) can be described as 


ey = c- Ey, 


(2) 


where S,, is the neutron binding energy of the compound 
nuclei, Hy; is ith cascading y-ray’s energy, and Ep is 
the kinetic energy of an incident neutron in the center 
of the mass system. In this work, Sn = 7.89 MeV and 
7.59 MeV for compound nuclei of Br and 8tBr respec- 
tively. However, only a single S» can be used in Eq. 2. 
Method from ref. [31] was adopted in this work, S» of the 
most abundant isotope in the sample ("°Br) was chosen 
and the abundance of the 81Br isotopes in the sample is 
scaled according to its S, value. The deviation caused 
by the chosen of S» is less than 3% in this case. The ma- 
nipulation mentioned above is implemented by so-called 
weighting function (WF) [30]. Detector count of each de- 
positional energy is multiplied by WF to correct for the 
original y-ray efficiency of the detectors, so that the re- 
lationship expressed in Eq. 1 is achieved[24, 30]. In this 
work, energy and resolution calibrations for each individ- 
ual CgD¢ detector were performed using standard y-ray 
sources, including ®°Co and !37Cs. A realistic Monte 
Carlo simulation of the experimental conditions and the 
target setup was performed using the Geant4[32] toolkit 
and was used to obtain accurate response function. Fi- 
nally, a least-squares method minimizing 


x? = 5 (ce, = 


i 


> R(Ea, Byi)WF(Ea)àEa) , (3) 
By 


was used to determine the WF, where i is the group of 
mono-energetic y-rays used in the simulation, Ez is the 
low deposited energy threshold of the CgD¢ detector, and 
R(Ea, Eyi) is the detector response function. The WF 
of the natural KBr and gold samples were determined 
independently using this method. A detailed study of 
the possible sources of systematic uncertainties revealed 
that the PHWT has an accuracy of 3%[33]. 
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TABLE 1. Characteristics of the samples and experimental duration of our experiments. 


Experimental eq-egduration (h) 


Target Mass (g) Diameter (mm) Thickness (mm) Purity (%) Without filter With Alter 
TST AU 4.83+0.01 40.00.02 0.20+0.02 299.99 11.0 4.0 
KBr 11.17+0.01 30.00.02 2.00.02 299.9 11.0 4.0 
Ph 13.93+0.01 40.0+0.02 0.98+0.02 299.9 8.0 = 
a 2.86+0.01 40.00.02 1.02+0.02 299.9 8.7 = 
empty = = = — 8.0 6.0 
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Fig. 4. (a) Original efficiency of the CgDe detector; (b) Effi- 
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ciency of the detector weighted by the WF; (c) Coefficient of 192 
the weighted efficiency vs incident y-ray’s energy; (d) WF of 193 


the CgDe detector system with KBr and Au sample. 


B. Double-bunch unfolding method 


At present, the accelerator complex of CSNS oper- 
ates in the normal mode, each pulse contains two proton 
bunches, the interval between the two bunches is 410 
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energy resolution. 

In the normal mode of CSNS, the two bunches in each 
single beam pulse are basically identical, and the tem- 
poral structure is extremely reproducible from pulse to 
pulse; the details can be found in Refs[17, 34]. For the 
statistical time count spectrum measured in the normal 
mode, the counts of each time bin E; measured in such a 
mode should theoretically depend on the counts C; mea- 
sured in the single-bunch mode, the transformation from 
the single-bunch mode spectrum to the normal mode 
spectrum can be written in the form of a matrix as 


E=R.-C, (4) 


where the vectors E = (Eo,--- , Ei,- , En) and C = 


(Cos ,Ci,--- ,Cy) as well as the transformation ma- 
trix R can be expressed as 
Ray 
R22 
R Ra+4)ı ` (5) 
R+4)2 Rii 
Ratai » Ran 


where A is the number of bins corresponding to the offset 
of the double-bunch interval (410 ns), Rij = 4 if i = j or 
j =i—A, otherwise, Ri; = 0. Using Eq. 5, the unfold- 
ing problem can be treated as an inverse matrix problem, 
however, the inverse matrix method occasionally yields 
several negative-value bins and oscillations caused by the 
statistical uncertainty in the measurements. Thus, an 
iterative algorithm using Bayesian estimation was devel- 
oped, the inverse matrix problem can be replaced by an 
iterative process[34, 35], that is 


Ch = Ft 
C 


where c™ indicates the result of the kth iteration, and 


ns|17]. The neutrons generated by two proton bunches 199 the initial co) is set as the E; measured in the double- 


overlap with each other. Thus, the neutron energy res- 
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bunch distribution. Part of unfolding results of KBr tar- 


olution of the TOF measurement is reduced, especially 201 get is illustrated in Fig. 5 
in the higher-neutron-energy region (above hundreds of 202 
eV). An unfolding method based on Bayes’ theorem was 203 and the double-bunch unfolding method were normalized 
170 developed by H. Yi et al.[34] to obtain better time and 203 by the proton beam number and are shown in Fig. 6. 
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Fig. 5. (a) Comparison between the weighted original spec- 
trum of the Br sample (black curve) and the weighted original 
spectrum obtained from the unfolding process (red curve); 
(b) ENDF/B-VIII.0 database of Br. In panel (a), the black 
curve shows the presence of clear double peaks caused by the 
double-bunch beam structure, while the resonance structures 
are restored when using the unfolding procedure, as indicated 
by the red curve. 
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Fig. 6. Preprocessed and normalized (according to the proton 
beam number) original spectra of ™**Br, "**C, "**Pb, and the 
empty target. 
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C. Background 


207 There are two kinds of components contributing to the 
20s background level in the capture cross section measure- 
209 ments with CgD¢ detectors[36]: sample-dependent back- 
210 ground Bsampie(tn) and sample-independent background 
211 Bempty (tn), that is 


B(ty) = Bempty (tn) =F Bsample(tn)- (7) 


213 The contribution of Bempty(tn) can be directly mea- 
214 sured with empty sample in same experimental condi- 
215 tions. On the other hand, the sample-dependent back- 
216 ground Bsample(tn) is caused by interactions between the 


212 


sample and all the types of in-beam particles, includ- 
ing the scattered-neutron-induced background By, (tn), 
scattered in-beam y-rays background Bs,(tn), and sam- 
ple activation background Bac.. Thus, the sample- 
dependent background can be expressed as[36] 


Bsample(tn) = Bsn(tn) 4 Bsy (tn) F Bag (8) 


The scattered-neutron-induced background Bsn(tn) can 
be determined by carbon sample measurement([39] 
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(W.Cc (tn) = W.Cempty(tn)) , (9) 


where Yc, e and Ypr, e1 are the neutron scattering yields 
of the carbon and Bromine targets obtained from the 
database, W.Cc(tn) and W.Cempty(tn) are normalized 
weighted counts of carbon sample and empty sample. 
The in-beam y-rays originated from neutron captures in 
the water moderator of the spallation source. Indeed, 
these y-rays can be scattered by the sample. The tar- 
get and energy dependence of in-beam y-ray background 
components were determined from a measurement of a 
lead sample and the absorption valleys of 4.28 eV, 5.18 
eV, 132 eV and 5.02 keV of the Ta—Ag—Co filter, as 
shown in Fig. 7. In this figure, the empty background 
Bempty(tn) is subtracted from all spectra, and the back- 
ground due to the scattered neutrons from lead sample 
is subtracted using Eq. 9. The figure also shows the ac- 
tivation background, which is determined by fitting the 
spectra’s platform above 11 ms (E, œ 0.2eV) [36]. In this 
region, neutrons are absorbed by the cadmium absorber, 
and the in-beam y-rays can be ignored; the counts in the 
residual TOF spectrum are attributed to the activation 
of the sample and the surrounding materials. While, in 
Fig. 7 (a), absorption valley of 5keV didn’t match the in- 
beam gamma-ray background shape determined by Pb 
target. We believe that is caused by delayed gamma-ray 
from inelastic scattering reaction channel of ®'Br. 81Br 
contains a level 536.2keV with lifetime of 34.6 us. In- 
elastic scattering reaction channel of ®'Br will open with 
neutron energy over 500keV ( 15 us of flight time for 
Back-n facility). Flight time of 5 keV absorption valley is 
about 80us. In region of 5 keV absorption valley, delayed 
gamma-ray form ®!Br will contribute some extra gamma- 
ray counts. While flight time for 132 eV absorption val- 
ley is 460 us, far more than lifetime of 536.2keV level. 
So, we use 132 eV black resonance to determine the in- 
beam gamma-ray background for KBr crystal. The level 
536.2keV will deexcite through gamma-rays of 260.2 keV 
and 276 keV. Thus, a threshold of 300 keV can be used 
in experimental data to verify this opinion, the results 
265 were shown in Fig. 7 (b). In this figure, absorption val- 
266 leys of both 132 eV and 5.02 keV agree with the in-beam 
267 gamma-ray background shape determined by Pb target. 
26s However, majority of gamma ray energy released by the 
269 composite nuclear decay generated by neutron capture 
270 is also below 300keV, so this threshold was not used in 
271 the experimental data processing process. 
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Fig. 7. (a) The measured in-beam y-ray backgrounds were 
normalized to match the values of the energies of the black 
resonances, namely 4.28 eV (Ta filter), 5.18 eV (Ag filter), 
and 132 eV and 5.016 keV (Co filter). (b) Used a threshold 
of 300 keV to avoid the influence of the 536 keV level. 


D. Experimental corrections and absolute neutron flux 
normalization 
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274 Inthe RRR, sample-related corrections were included 
275 in the SAMMY(26] analysis. In the URR, neutron mul- 
276 tiple scattering events and self-shielding corrections in 
277 the sample were determined through Geant4 simulations, 
278 that is 


Nw Ot 


ISa ` 1 — exp (N feto)” (10) 


279 Oy = 


230 where fe is the correction factor for the multiple scat- 
281 tering events and the self-shielding effects, p is the 
232 area density in atom/barn of the sample, which is 
283 0.00339 atom/barn for the Bromine sample and 0.00117 
234 atom/barn for the gold sample. Uncertainty of fe fac- 
ess tor be considered as 1%. A relative normalization of the 


287 well-defined energy dependence of the neutron flux could 
288 be obtained from the various runs via the Li-Si neutron 
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Fig. 8. Correction factor for the multiple scattering events 
and the self-shielding calculated using the Geant4 toolkit: 
the black solid dots represent natural Bromine, while the red 
solid dots represent the gold sample. 
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Fig. 9. Capture yield of the first resonance of 1°’ Au (4.9 
eV) measured at the Back-n facility and normalized by the 
SAMMY[26] fit. 


detectors in ES#1. The absolute flux was determined by 
means of the gold reference sample using the (n,y) cross 
section of 197 Au as a standard. 


The first gold resonance at 4.9 eV was used to de- 
fine the flux in the RRR using the saturated resonance 
method, as shown in Fig. 9. The absolute yield nor- 
malization was determined through a fit of gold target’s 
data using the R-matrix code SAMMY [26] and by adopt- 
ing the resonance parameters of Ref. [39]. A systematic 
uncertainty of 1.5% has been adopted for the absolute 
299 flux normalization. In the keV region, the average (n,y) 
300 cross sections were obtained relative to gold. The back- 
301 ground of the Au spectrum was determined using the 
302 Same method adopted for the Bromine spectra. 
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E. Discussion of the uncertainties 


The total uncertainties, including statistical and sys- ‘ 


tematic, will be discussed. The statistical uncertainty 
comes from raw counts in a energy bin of four samples 
and was estimated to be <2.0%. In fact, since the raw 
counts will change depending on the width of energy bins 
and value of (n, y) cross sections. Wider energy bins will 
help to increase the counts and reduce the statistical er- 
ror, for energy >2.0 keV), but at the same time, too wide 
energy bins cannot exhibit the fine resonances structure. 

The systematic uncertainty are mainly contributed 
by the uncertainty of experimental conditions and data 
analysis method. Uncertainty of experimental conditions 


contain several types including uncertainty of sample pa- ° 
rameter, beam profile on sample, neutron energy spec- : 
trum and proton beam power. Uncertainty of data anal- : 


ysis method are mainly caused by the PHWT method, 
double-bunch unfolding process, normalization, back- 
ground subtraction and correction in the URR. Finally, 
according to the error propagation, the overall experi- 
mental uncertainty is less than 10.60% in Table 2. Such 
a high error mainly comes from the uncertainty of the 
neutron spectrum (<8%). Therefore, a good neutron 
energy spectrum with lower uncertainty will greatly im- 
prove the accuracy of this experiment. 


IV. RESULTS AND DISCUSSION 
A. R-matrix fits 


In the region from 1 to 2000 eV, these capture 
yields were analyzed using the R-matrix analysis code 
SAMMY 26]. The yield is parametrized via Reich-Moore 
approximation to R-matrix formalism. Scattering radius 
of 6.85 fm and a temperature of 293 K for the correction 
of the Doppler effect were adopted. Other experimental 
effects, that is, multiple neutron scattering in the sam- 
ple and neutron self-shielding, are properly taken into 
account within the SAMMY code. The resonance broad- 
ening due to the neutron energy resolution function was 


mented RPI parametrization|[40]. 

The fitting procedure allowed us to extract the res- 
onance parameters (radiation width Ty, neutron width 
Tn, orbital angular momentum L, etc.) from the mea- 
sured capture yields. However, in many cases, only 


the resonance energy Fr and the total capture kernel ie 


k should be considered as the real measurable quanti- 


ties. The capture kernel k is proportional to the area 3 


under an isolated resonance and is given by 


k = gVpVy/T, (11) 


where g is a statistical factor defined as 


2J +1 


= E E = 


also considered in the SAMMY fits through the imple- | 
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399 


400 


40 


pa 


where J is the total angular momentum, s is the spin of 
the incident particles (s = 4 for neutrons), and T is the 
spin of the target particles. 

Resonance parameters of Bromine contained in 
ENDF/B-VIII database were adopted for initial param- 
eters of SAMMY fitting procedure. Resonance param- 
eters from ENDF/B-VIII database is consistent with 
measurement of R. L. Macklin (1988)[14], but con- 
tains more small resonances than result of M. Ohkubo 
et al(1981)[13]. The final SAMMY fitted results of the 
KBr crystal capture yield are shown in Fig. 10: the 
black data represent the experimental capture yield mea- 
sured in this work, and the red solid curve is the actual 
SAMMY fit to the present data. The fitted resonance pa- 
rameters and the radiative kernels derived using Eq. 11 
are listed in Table 5. No external R-function was em- 
ployed. The contributions of 394K were considered with 
resonance parameters from ENDF/B-VIII database. 

For comparison with evaluated databases, we calcu- 
lated neutron capture cross section Jexp.(n, Y) of Br with 
our experimental yield Yoxp.(n, Y), 


Yoxp. (0, Y)Ctot 
(1 = e— Norot) 


in which Gtot is the total cross section calculated from 
the Evaluated Nuclear Data Files (ENDF/B-VIII.0) 
database, and N is atomic number density of Br in KBr 
crystal. Calculated oexp.(n, Y) is plotted with blue dots 
in Fig. 12. 


Cexp. (HV) = (13) 


B. Statistical analysis 


The present set of resonance parameters was used for 
a statistical analysis to determine the nuclear proper- 
ties required for the cross-section calculations[27]. The 
cumulative number of resonances as a function of neu- 
tron energy is shown in Fig. 13 (a) and (b) for Br and 
81Br respectively. This figure provides an efficient way 
to investigate level populations and missing levels. The 
average s-wave level spacings Do are directly related to 
the inverse slope of these plots and can be obtained by 
the linear least-squares fits indicated by straight lines, 
which is 57.397 eV and 27.855 eV for Br and ®!Br re- 
spectively. Points fall below the fitted straight line indi- 
cating that levels are being missed in the resonance anal- 
ysis. In addition, the average radiative widths < Ty > 
are calculated with SAMMY fitted resonance parameters 
which is 287.3410.3 meV and 295.5+12.5 meV for Br 
and °'Br respectively. 


C. MACSs 


In continuum region below 300 keV, average cross 
sections were obtained with resolution of 20 bins per 


Component 


Uncertainty (%) 


PHWT 

Unfolding method 
Normalization 

Background subtraction 
Experimental corrections 
Proton beam power 

Neutron beam profile 

Target parameters 

Neutron spectrum (> 0.15 MeV) 
Neutron spectrum (< 0.15 MeV) 
Statistical 


3.0 
3.0 
1.5 
2.0 
1.0 
1.5 
1.5 
0.1 
4.5 
8.0 
2.0 


Total uncertainty(> 0.15 MeV) 
Total uncertainty(< 0.15 MeV) 


8.2 
10.6 


TABLE 2. Statistical and systematic uncertainties of the experiment. 
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Fig. 10. Analysis of the resonance parameters of the experimental data fitted by the R-matrix code SAMMY. 
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function of the resonance energy. (b) Corresponding distribution of the ratios. 


(a) Ratio of the capture kernels obtained from the ENDF/B-VIII.0 database and the present C6De results as a 
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Fig. 12. Capture cross section exp. (n, y) calculated from the experimental yields obtained in this work (black dots), SAMMY 
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Fig. 13. Staircase plots of the cumulative numbers of resonances in the investigated Bromine isotopes. 


402 decade. Averaged cross section relative to the gold sam- 
a03 ple ogr(En) is given, according to 


Caul En) 
(Tau (En)}) 


where (oau (En)) and (opgr(En)) are the experimental val- 
ues measured in this work, and oau(Ey) is the evaluated 
value obtained from the ENDF/B-VIII.0 database. The 
uncertainty of the standard 1% Au cross section recom- 
mended by the ENDF/B-VIII.0 library was estimated to 


404 Opr(En) = (oBr(En)) ; (14) 


aio be in the order of 6.0% below 200 keV and 4.0% between 
41 200 and 300 keV. In URR, average cross section of K 
a12 is lower than Br with 2 orders of magnitude, neutron 
a3 Capture events of K atoms can be ignored in this region. 
414 Finally, the results of the cross section of Br in the con- 
ais tinuum region from 10 to 300 keV are given in 20 bins 
ais per decade and are listed in Table 3. 

In Fig. 14, the averaged cross sections obtained in 
ais this work in the continuum region(red dots) are com- 
aio pared with previous experimental results and the eval- 


417 


420 


421 


422 
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TABLE 3. Averaged capture cross sections of natural Bromine in the URR, op; calculated with Li-Si neutron spectrum also 


gived in this table. 


Fiow (keV) Eup (keV) opr (Au) (mb) opr (Li-Si) (mb) 


Eiow (keV) Eup (keV) opr (Au) (mb) op, (Li-Si) (mb) 
2.0 23 2554.0£247.6 2554.0£247.6 
2.3 2.6 3063.6+291.7 3063.6+291.7 
2.6 2.9 2642.6+250.4 2642.6+250.4 
2.9 3.3 2475.0+234.6 2475.0+234.6 
3.3 3.7 1834.4+173.7 1834.4+173.7 
3.7 4.2 2146.6+203.4 2146.6+203.4 
4.2 4.8 2261.1+214.2 2261.14+214.2 
4.8 5.4 1251.1+118.6 1081.1-+118.6 
5.4 6.1 1314.9+124.4 1511.4+160.1 
6.1 6.9 1445.7+136.6 1389.3+139.6 
6.9 7.8 1516.2+143.2 1543.1-+157.9 
7.8 8.9 1249.5+117.9 1282.5+131.5 
8.9 10.1 1300.4+122.5 1416.5+148.9 
10.1 11.4 1141.4+122.5 1175.2+121.3 
11.4 12.9 1211.3+107.3 1223.5+127.5 
12.9 14.6 874.8+78.8 824.1+88.9 
14.6 16.5 875.2+82.1 882.8+92.4 
16.5 18.7 840.1+78.5 847.1+89.1 
18.7 21.2 845.5+78.8 863.8+91.1 
21.2 24.0 738.6£68.8 714.2+74.5 
24.0 27.2 575.2+68.8 A75.54£52.7 


27.2 30.8 658.5+61.1 728.3£77.3 
30.8 34.9 527.1448.7 477.0+51.5 
34.9 39.5 560.1-£51.7 588.4+61.2 
39.5 44.7 496.1+45.6 490.353.3 
44.7 50.6 430.3439.5 398.6£41.5 
50.6 57.3 438.8+40.2 435.6+47.5 
57.3 64.9 368.7+433.6 332.6436.9 
64.9 73.5 337.6£30.7 316.0+34.4 
73.5 83.3 310.7428.2 283.7+30.2 
83.3 94.3 335.6£30.3 348.1+38.2 
94.3 106.8 302.927.3 308.531.3 
106.8 120.9 273.8+24.6 270.5+29.7 
120.9 136.9 236.4+21.1 213.4+24.1 
136.9 155.0 251.8+22.4 251.8+28.3 
155.0 175.5 243.0+21.5 251.9+27.9 
175.5 198.7 209.3+18.4 191.5+20.9 
198.7 225.0 204.1+17.9 202.9+21.1 
225.0 254.8 186.0+16.2 200.7+21.0 
254.8 288.5 152.5+13.2 147.9+15.1 
288.5 326.7 131.7+11.4 154.9+15.3 
326.7 370.0 112.8+9.7 133.5414.2 
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Fig. 14. Capture cross sections of Bromine in this work rel- 447 
ative to the 1°” Au sample (red square dots). Cross section 448 
calculated with neutron flux determined by Li-Si detector is 
also plotted for comparison. 


uated database. Our measurements are consistent with 452 
the results of H.Gibbons et al. (1961), while higher than 453 
all evaluated database and results of Ju. P. Popov et al. 454 
(1961). For comparison, in this figure we also plotted 455 
(black dots). In region over 200 keV, Li-Si detector de- 456 


average cross sections in the URR. The calculations 
were based on the Hauser—Feshbach statistical emission 
model, which assumes that the capture reactions occur 
by means of a compound nuclear system that reaches a 
statistical equilibrium. The previous obtained statistical 
average level space Do average radiation width < I, > 
were used as input parameters for the TALYS code cal- 
culations. In addition, the global neutron optical model 
potential of Ref. [41] was used in the calculations and 
other parameters are chosen with method reported in 
S. Chen et al. (2020)/[38], photon strength function is 
given by Kopecky and Uhl[43, 44], level density a and 
nuclear temperature T are given by Gilbert-Cameron 
model with adjusted parameters. In addition, < Ty > 
is systematically multiplied by a factor of 0.9 for both 
isotopes to obtain y transmission coefficient. The calcu- 
lated capture cross sections well reproduced the available 
experimental average cross sections of Br and ®!Br as 
illustrated in Fig. 15 (a) and (b). 


Fig. 16 shows the comparison of our results, the evalu- 
ated databases, and the recommended values compiled in 
the Karlsruhe Astrophysical Database of Nucleosynthe- 
sis in Stars (KADoNiS): (a) the MACS values of ?Br ob- 
tained in this work, which are basically located between 
the values of the JEFF-3.3, TENDL-2021, ENDF/B- 
VIII.0, and JENDL-5 databases, are in good agreement 
with the KADoNiS values; (b) for ®!Br, the calculated 


a25 termined result is obviously higher than result measured 457 values are considerably higher than those of the evalu- 


426 


427 


428 


429 


relative to gold sample. In addition, fluctuations near 30 45s 
keV of results determined by Li-Si detector is caused by 
aluminum material in neutron beam pipe. 465 


The TALYS code was used to describe the isotopic 460 


ated database and the KADoNiS recommended values. 


From these averaged cross section values, the MACSs 
of “Br and 81Br were calculated for thermal energies kT 
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Fig. 15. Averaged capture cross sections of Br and ®!Br calculated using the TALYS code. The theoretical averaged cross 
sections of “Br are basically consistent with the previous measurements and evaluated values, while those of *'Br roughly 
agree with the previous empirical values but are higher than those of the evaluated database. 


TABLE 4. MACSs of Br and ŝ'Rb. 


kT (keV) Br (mb) 8TBr (mb) 
This work(mb) KADOoNiS v1.0 This work KADOoNiS v1.0 
5 1938+193 1890+181 725+72 71535 
10 1272+127 1223+105 53453 44727 
15 1003100 96674 433443 357+21 
20 852+85 823+58 369+36 307+16 
25 15375 729+49 32532 272+13 
30 682+68 661+44 293+29 248+10 
40 585+58 567+38 248+24 212+10 
50 519+51 50335 21821 188+10 
60 47047 45433 196+19 170+10 
80 396+39 383+30 165+16 145+9 
100 34134 332+27 143+14 127+8 
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Fig. 16. MACSs calculated with Eq. 15 are illustrated in (a) and (b) for ”°Br and ®'Br, respectively. The blue triangles are 
obtained from ref. [3] 


461 


from 5 to 100 keV according to 


2 f o(En) Ene TAE, 
PES [5° Ene~2o/ FTE, i 


(15) 


and the corresponding results are listed in Table 4 (a) 
and (b), respectively. Fig. 16 shows the comparison of 


s our results, the evaluated databases. 


Fig. 16 shows the comparison of our results, the eval- 
uated databases, and the recommended values compiled 
in the Karlsruhe Astrophysical Database of Nucleosyn- 
thesis in Stars (KADoNiS): (a) the MACS values of 7°Br 
obtained in this work, which are basically consistent with 
the KADoNiS values and evaluated databases; (b) for 
81Br, the calculated values are considerably higher than 
those of the evaluated database and the KADoNiS v0.3 
and v1.0 recommended values but smaller than KADo- 


s NiS v0.0. 


Fig. 17 shows the comparison of our results with the 
previously recommended MACSs at the nuclear astro- 
physics concerned thermal energy of kT = 30 keV. The 
value of 638464 mb for Br showed in Fig. 17 (a) is in 
good agreement with the KADoNiS recommended value 
of 661+44 mb within the uncertainty range. On the 
other hand, the MACS value of 293429 mb for 8tBr 
shows a clear discrepancy with the KADoNiS recom- 
mend value of 248+10 mb. 


V. CONCLUSIONS 


(n,y) reaction of natural Bromine were measured at 
the Back-n facility using an array of four CgDg¢ detec- 
tors. The PHWT with Monte-Carlo simulations and the 
double-bunch unfolding method were used for the data 
preprocessing. Black resonance method with a Ta—Ag— 
Co filter and dedicated measurements were used to study 
the experimental backgrounds to obtain accurate back- 
grounds. 

The capture yields have been analyzed in resolved res- 


s onance region with the R-matrix code SAMMY. A to- 


tal of 121 resonances were found in the neutron energy 
range from 1 to around 2000 eV. From these results aver- 
age level spacing, radiative widths, and neutron strength 
functions have been deduced by statistical analyzes to es- 
tablish a consistent set of input data for detailed cross 
section calculations with the Hauser-Feshbach statistical 
model. The MACSs for 7Br obtained in this work are 
basically located between the JEFF-3.3, TENDL 2021, 
and ENDF/B-VIII.0 JENDL-5 databases, and they are 


s in good agreement with the KADoNiS v1.0 values. By 


contrast, for 8tBr, the calculated values are substantially 
higher than those of the evaluated database and the 
KADoNiS recommended values. The MACSs at kT = 30 
keV are 638464 and 293429 mb for “Br and ®!Br, re- 
spectively. Our results provide extra constraints on the 
actual MACSs of ’’Br and ®!Br. 
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Fig. 17. MACSs at the thermal energy of kT = 30 keV for (a) “Br and (b) "Br. The results from the present work are indicated 
by the yellow-shaded areas. The black dot represent previous experimental results: HKUO08b[3], Mac88a, WBK86b/42], 
AGM71[45], HeP06[2], Gor05[47], RaT99[48], Har81[46], WFH78([16] and the dots labelled endfb8, jeff33, tendl5, jendl21 are 
calculated from the evaluated databases according to Eq. 15. The red-triangle dots labelled AGM71 represents the evaluated 
value taken from Ref. [37], while the three dots labelled KADoNiS v0.0, KADoNiS v0.3 and KADoNiS v1.0 are recommended 
values from three version of KADoNiS database. 


TABLE 5: Resonance parameters extracted via SAMMY/([26] fit to our experimental data. 
ENDF/B-VIIL.0[8] are still listed for comparison. 
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Resonance parameters from 


Ga this work ENDF/B-VIIL.0[8] aes J this work ENDF/B-VIIL.0[8} 
nu al 

Eexp. [eV] Kexp. [meV] EenatleV] kenat[meV] Eexp. eV] kexp. [meV] EenatleV] Kenar[meV] 
Br 2 35.82+0.015 12.08-40.52 35.8 22.5 51Br 2 1106.0340.246 98.75+3.49 1103.0 103.44 
Br 1 53.74+0.018  10.67-40.62 53.7 10.9 Br 2 1140.99+0.439  3.58+0.68 1138.0 3.4 
SBr 2 101.250.039 44.08+1.66 101.2 56.83 S1Br 1 1147.39+0.236 93.47+4.6 1147.0 89.5 
S1Br 1 135.64+0.053 50.08-+6.31 135.6 65.5 "Br 1 1165.64+0.027 4.83+0.85 1165.0 4.8 
Br 2 158.93+0.063 0.24+0.05 157.9 0.3 Br 1 1186.49+0.278 15.89+2.79 1187.0 20.5 
Br 1 189.59+0.037 21.4+0.94 189.5 23.9 Br 1 1192.33+1.153 1.17+0.23 1192.0 1.2 
Br 2 192.64+0.073  1.11+0.11 192.5 1.9 Br 2 1201.44+0.33 164.81+12.11 1201.0 177.9 
S1Br 1 205.17+0.037 5.94+0.26 205.0 5.41 S1Br 1 1207.23+40.806 88.0-+13.31 1209.0 84.55 
Br 2 211.78+0.069 0.60.04 211.6 0.6 Br 2 1228.64+0.233 29.02+1.9 1228.0 21.2 
Br 2 238.7640.051 117.32+4.64 238.9 120.9 Br 2 1236.61+0.141 1.97+0.39 1239.0 2.0 
S1Br 2 255.95+0.102  0.46+0.04 255.6 0.55 51Br 2 1276.64+0.364 85.82+3.43 1276.0 147.95 
Br 1 292.13+0.471 0.65+0.12 292.5 0.6 Br 2 1293.52+0.093  4.39+0.86 1296.0 4.4 
Br 1 294.480.058 23.6841.77 294.3 22.1 Br 2 1300.82+1.237 0.74+0.15 1301.0 0.8 
Br 2 318.84+0.074 124.39+5.36 318.6 122.7 Br 2 1312.89+0.182 27.214+2.86 1312.0 23.7 
SIBr 2 336.8740.207 0.78+0.1 336.7 0.74 ®1Br 3 1312.1741.176 2.31+0.44 1312.0 2.25 
S1Br 1 341.07+0.543 0.13+0.02 340.9 0.13 "Br 2 1317.76+0.291 2.40.47 1317.0 2.4 
S1Br 2 347.92+0.452  0.29+0.06 348.2 0.3 S1Br 3 1340.78+0.057 2.68+0.53 1342.0 2.7 
S1Br 3 369.44+0.108  1.39+0.1 369.3 1.13 Br 2 1349.2+1.371 0.49+0.1 1349.0 0.5 
Br 2 395.310.066  36.84+3.29 394.6 47.4 Br 2 1358.55+0.092 1.68+0.33 1362.0 1.7 
Br 2 465.24+0.243 1.06+0.21 464.2 1.1 Br 1 1379.84+0.177 59.27+4.18 1380.0 31.01 
Br 1 468.79+0.103 20.44-+1.52 468.2 26.2 Br 2 1390.9+1.417 2.16+0.38 1391.0 1.9 
Br 1 482.98+0.132  22.0+1.81 482.7 27.7 "Br 2 1416.77+1.096 4.71+0.69 1415.0 3.3 
Br 2 491.55+0.25 0.94+0.13 490.8 0.9 S1Br 3 1442.5+0.07 3.88+0.76 1441.0 3.86 
Br 2 510.28+0.525 0.25+0.05 510.2 0.3 Br 2 1447.66+1.488 2.45+0.48 1448.0 2.4 
Br 1 548.78+0.226 1.56+0.16 548.2 1.3 Br 2 1454.83+0.268 159.58+8.84 1455.0 128.3 
SIBr 1 560.05-40.285  6.95-+0.99 560.2 6.89 Br 2 1463.9141.425 9.92+1.83 1464.0 9.6 
Br 2 564.98+0.066 159.43+9.44 564.9 117.7 Br 1 1469.54+0.236 126.27+8.53 1470.0 81.7 
®1Br 2 578.62+0.319 10.98+1.48 578.7 94.1 "Br 2 1483.27+0.337 22.05+1.51 1483.0 10.8 
Br 1 604.79+0.069  74.95+2.2 604.0 78.0 Br 2 1531.28+0.482 110.77411.92 1531.0 170.7 
Br 1 638.07+0.083 42.04+3.25 637.9 33.2 ®1Br 2 1543.09+0.484 164.81+12.06 1548.0 128.81 
Br 2 645.92+0.076  104.0+3.18 646.2 98.83 "Br 2 1572.06+0.332 67.9345.44 1572.0 53.2 
S1Br 2 668.57+0.097 94.42+2.75 668.5 133.86 | Br 2 1590.14+0.211 168.89+9.75 1590.0 151.1 
Br 0 707.04+0.146 0.58+0.11 708.0 0.59 "Br 1 1630.24+0.098  6.68+1.19 1634.0 5.81 
Br 2 749.73+0.084 112.09+9.76 749.7 89.2 Br 2 1648.59+0.113 2.54+0.5 1651.0 2.5 
SIBr 1 771.87+0.259 0.86+0.12 771.8 0.67 S1Br 2 1667.03+0.312 2.69+0.53 1666.0 2.74 
Br 2 789.0340.1 101.74£16.57 788.3 138.1 Br 2 1671.6440.089 3.95+0.77 1674.0 3.9 
Br 1 800.2+0.674  0.69+0.15 800.7 1.2 Br 1 1688.9+1.575 2.26+0.39 1686.0 2.0 
Br 2 815.0140.594 1.67+0.25 814.1 1.5 S1Br 1 1707.18+0.676 34.7344.43 1708.0 35.63 
Br 1 819.85+0.678  1.06+0.18 818.4 0.9 Br 1 1718.49+0.534 112.25+13.78 1720.0 87.1 
Br 2 832.52+0.095 21.3641.59 831.7 23.2 "Br 1 1720.8+0.581 75.28+9.91 1723.0 65.9 
S1Br 1 850.9+0.1 20.86+2.01 850.2 20.11 Br 1 1734.14+0.951 1.60.32 1734.0 1.6 
Br 2 871.5+0.089  4.28+0.83 870.2 4.9 Br 2 1746.83+0.064 10.8+2.01 1744.0 10.4 
Br 2 893.16+0.11 23.97-+1.25 892.7 21.8 Br 1 1772.2+0.328  40.79+2.1 1772.0 86.9 
Br 2 915.83+0.14 1.31+0.26 914.8 1.3 Br 2 1782.61+1.36  3.66+40.74 1785.0 4.3 
Br 1 931.6640.108  92.9+7.94 930.5 72.5 Br 1 1797.380.078 8.95+1.68 1803.0 9.2 
8&iBr 1 961.980.079 1.23+0.23 959.9 1.2 S1Br 2 1822.540.284 10.83+41.92 1824.0 10.84 
SIBr 1 994.59+0.125 33.02+2.39 994.0 20.65 Br 2 1829.69+1.0 165.65+25.0 1829.0 172.5 
Br 2 1009.96+0.722 1.78+0.29 1012.0 1.5 51Br 2 1834.21+0.571 86.82+14.21 1834.0 96.19 


continued on next page 
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(continued TABLE.5) 


coed this work ENDF/B-VIIL.0[8] T this work ENDF/B-VIIL.0[8] 
Exp. [eV] Kexp.[meV]  Eenat[eV] Kenat[meV] Exp. [eV] kexp.[meV]  Eenar[eV] kena¢[meV] 
Br 1 1025.15+0.295 7.8840.59 1025.0 5.6 Br 2 1872.34+0.573 86.24+4.93 1874.0 191.8 
Br 2 1029.89+0.191  1.2+0.24 1031.0 1.2 ®1Br 1 1891.79+0.058 30.41+4.45 1897.0 29.44 
SIBr 3 1039.82+0.417  2.81+0.55 1038.0 2.87 S'Br 1 1904.58+0.067 3.67+0.65 1907.0 3.66 
Br 2 1043.6540.158 29.0641.66 1043.0 25.7 Br 2 1942.35+0.073 4.62+0.9 1948.0 4.6 
SIBr 3 1069.040.221 7.85+0.41 1069.0 4.51 SIBr 3 1953.94+0.432  2.52+40.5 1952.0 2.53 
SIBr 1 1083.89+0.142 13.0741.98 1082.0 14.57 Br 2 1969.1740.256 161.57413.28 1969.0 108.2 
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